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ABSTRACT 

Aims. The aim of the work presented in this paper is to test and optimise supernova detection methods based on the optimal image 
subtraction technique. The main focus is on applying the detection methods to wide field supernova imaging surveys and in particular 
to the Stockholm VIMOS Supernova Survey (SVISS). 

Methods. We have constructed a supernova detection pipeline for imaging surveys. The core of the pipeline is image subtraction using 
the ISIS 2.2 package. Using real data from the SVISS we simulate supernovae in the images, both inside and outside galaxies. The 
detection pipeline is then run on the simulated frames and the effects of image quality and subtraction parameters on the detection 
efficiency and photometric accuracy are studied. 

Results. The pipeline allows efficient detection of faint supernovae in the deep imaging data. It also allows controlling and correcting 
for possible systematic effects in the SN detection and photometry. We find such a systematic effect in the form of a small systematic 
flux offset remaining at the positions of galaxies in the subtracted frames. This offset will not only affect the photometric accuracy of 
the survey, but also the detection efficiencies. 

Conclusions. Our study has shown that ISIS 2.2 works well for the SVISS data. We have found that the detection efficiency and 
photometric accuracy of the survey are affected by the stamp selection for the image subtraction and by host galaxy brightness. With 
our tools the subtraction results can be further optimised, any systematic effects can be controlled and photometric errors estimated, 
which is very important for the SVISS, as well as for future SN searches based on large imaging surveys such as Pan-STARRS and 
LSST. 

Key words. Supernovae: general ~ Methods: data analysis ~ Surveys 



1. Introduction 

Imaging surveys to find and characterise variable sources are 
used to study a number of different astronomical phenomena. 
Detection and photometry of variable sources can be done by 
subtracting images taken at different epochs from each other 
and analysing the subtracted frames, a technique known as dif- 
ference imaging. In addition to supernovae (SNe), this tech- 
nique has been used by various authors investigating many dif- 
fere nt subjects, am ong which are variable stars in crowde d fields 
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Corwin et alj| 2006), gravitational micro-lens ing (e.g.lAlard 
Sumi et alj|2006.) and planet detection (e.g. lHolman et al 



20071) . 

The most complicated step in the difference imaging pro- 
cess is taking into account the spatial and temporal variability of 
the point spread functions (PSFs) of the images. The difference 
imaging technique used in this work is based on the Optimal 
Image Subtraction (OIS) code first presented in Alard&^Lupton 
d 19981 hereafter AL98), later improved upon in lAlardI (l2000l 
hereafter AOO) and available in the ISIS 2.2 package. In these 



papers the code is tested on variable stars (i.e. point sources with 
varying brightness with little or no background light). When de- 
tecting SNe the situation can be quite different. In this type of 
data the background light from the host galaxy often dominates 
the total flux at the location of the supernova. 

Difference imaging, and OIS in particular, has been used to 
detect SNe in many projects. Recently a number of large SN 
surveys, mainly aiming at finding thermonuclear SNe for deter- 
mining cosmological parameters, have released their first results. 
Image subtraction using the OIS algorithm followed by detection 
of SN candidates in the subtracted frames is a technique that has 
been used by the ESSENCE surv ey (the obser vations and data 
analysis are des cribed i n Miknaitis et al. '2007') , the Supernov a 
Legacy Survey ( Astier et al..,2006,) and Poznanski et alT(l2007l) . 
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ICappeUaroetal.1 ( l2005h and iBotticella et alj (l2008l) used flie 
same technique in their SN surveys aiming to determine the 
core collapse SN rate at intermediate redshift. Difference imag- 
ing h as also successfully been used in near-infrared SN searches 
(e.g. iMattila&Melki^boOU iMaiolino etaT]|2002h concentrat- 
ing on individual starburst galaxies. For surveys with the aim 
of finding supernova rates the detection efficiency is the most 
important issue, while the SN photometry is of somewhat less 
importance. However, for surveys targeting la SNe the photom- 
etry is of utmost importance since the photometric accuracy di- 
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rectly determines how well the cosmological parameters can be 
constrained. Therefore, in these programmes the photometry is 
treated very carefully and not necessarily obtained from the sub- 
tracted frames. For studies of individual SNe difference imaging 
has been used to obtain accurate photometry of supernovae with 
a complicated background. This is particularly important at late 
times when the SN has become too faint compared to the back- 
ground host galaxy in the unsu btracted image for PSF photome- 
try to give reliable results (e.g. ISoUerman et ani2002h . More re- 
cently, difference imaging has also successfully been applied to 
mid-infrared SN o bservations with the Spitzer Space Telescope 
jMeikle et al.1 '2006) and to adaptive optics assisted detection of 
SNe in luminous infrared galaxies (Mattila et al. 2007). 

In this paper we investigate several aspects of difference 
imaging with ISIS 2.2 when used to detect and study supernovae. 
Previous papers have akeady reporte d some test s done on the 
image subtraction (e.g., AL98; AOO; Israel et al. 2007) and on 
the SN detection process (e.g. Mattila & Meikle 2001.) . The sur- 
veys for thermonuclear SNe all have extensive sections where 
the methods to ascertain the photometric accuracy ( both in terms 
of st atistical and systematic errors) is described (e.g.l Astier et al.l 
12006 ). However, no systematic tests on the subtraction method 
applied to supernova data have been reported. The main purpose 
of this paper is to describe our supernova detection pipeline and 
to study how it performs on wide field imaging surveys, in par- 
ticular the Stockholm VIMOS Supernova Survey (SVISS). The 
SVISS is an imaging survey aiming t o find core-collapse super - 
novae at 0.1 < z < 1.0 (for details see Dahlen & FranssonI 19991) . 
The main scientific goal for this project is to find reliable core 
collapse supernova rate densities for this redshift range. We thus 
need to have detailed knowledge of the supernova detection ef- 
ficiency. In this project the photometric accuracy is important to 
correctly type the SNe based on colour and light curve informa- 
tion. Another goal of this paper is to investigate how supernova 
detection and photometry (with ISIS used as the image subtrac- 
tion tool) are affected by the selection of subtraction and detec- 
tion parameters and the presence of a host galaxy. In this paper 
we discuss the subtraction and detection method, in the follow- 
ing papers we will present the light curves of the detected SNe, 
the typing methodology and calculations of SN rates at high red- 
shifts. 

The first part of the paper (section |2]i contains a detailed de- 
scription of the standard image subtraction method used in this 
work. In section|3]we describe the data sets, in section|4]how the 
simulations of supernovae have been setup and in section |5] we 
present the results from the extensive testing. Finally, we discuss 
the results and possible implications in section |6] and conclude 
by giving a short summary in section Q The Vega magnitude 
system has been used throughout the paper. 

2. Image subtraction and SN detection method 

We have constructed a supernova detection pipeline suitable for 
large images and multi-epoch data. The pipeline consists of a 
number of IRAF/PYRAF scripts that are run in sequence, and go 
through the following steps: (i) accurate image alignment over 
the entire frame; (ii) convolving the better seeing image to the 
same PSF size and shape as the poorer seeing image, using a spa- 
tially varying kernel; (iii) subtracting the images; (iv) detection 
of sources in the subtracted frames, using both source detection 
software and eye-ball detection; (v) photometry and construction 
of light curves of the detected objects. 

The input to the pipeline consists of one reference image and 
a list of later epoch images. The reference image is subtracted 



from each of the later epoch images and objects with varying 
brightness are detected in the first search epoch image. The out- 
put from the pipeline is a list of supernova candidates detected 
in the first search epoch and a light curve constructed from mea- 
surements on the detected position in the subsequent epoch im- 
ages. The detection part of the pipeline can then be rerun with 
the previous second epoch as detection epoch (enabling SNe that 
have appeared between the first two epochs to be found). This 
can be repeated for all search epochs, i.e. once a particular search 
epoch has been run through the pipeline, that epoch can then be 
used as a reference epoch for subsequent epochs. This is use- 
ful when a search epoch image is of better quality (in terms of 
seeing and depth) than the reference epoch image. 

2. 1 . Registering images 

The subtraction method is quite sensitive to how well the im- 
ages are registered. Spatial variability of the convolution kernel 
can somewhat co mpensate for a no n-perfect image alignment 
(also discussed in llsrael et al.ll2007l) . but in general the images 
should be aligned with high precision before running the actual 
subtraction. Accurate image alignment is particularly important 
when the seeing differences are small, since narrow convolution 
kernels will not be able to compensate for poor image alignment. 

All the input images are aligned to the reference image by 
using the IRAF tasks geomap and geotran. A number of bright 
sources visible in all of the frames are used as input coordinate 
lists for these tasks. The exact positions of the sources are found 
by using the IRAF task imcentroid. The number of reference 
sources used is dependent on the order of the geometrical trans- 
formation and the size of the image. 

2.2. PSF matciiing and image subtraction 

The atmospheric conditions will affect the width and shape of 
the PSFs of the images. The combination of the atmospheric ef- 
fects with the optics, focusing and camera setup results in vari- 
ations of the PSF of the images both in spatial coordinates and 
over time (i.e. different epochs can have different PSFs, and also 
different spatial variation). A successful image subtraction tech- 
nique must be able to match the PSF of the reference image to 
the PSF of the later epoch images also allowing the convolution 
kernel to vary spatially. The OIS technique presented in AOO al- 
lows this. 

Following AOO (using the same notation) we summaries the 
OIS technique. The principal difficulty with image subtraction 
is to find a convolution kernel (K) that can transform the PSF 
of a reference image (R) to the PSF of an arbitrary image (/). 
The best-fit kernel can be found by minimising the sum over all 
pixels 

Y^i[R ® K](xi,yi) - I(xi,yd + bg(xi,yi))\ (1) 

where bg(xi,yi) is a spatially varying background term. The ker- 
nel can be written as a sum of basis functions 

K(u, v) = ^ a„(x, y)K„{u, v) (2) 

n 

where a„ contain the spatial variations of the kernel, which are 
polynomial functions of a given degree. K„ are the basis func- 
tions for the kernel and u and v denote the PSF kernel coordi- 
nates. The basis functions are Gaussian functions of the type 

Kn{u,v) ^ e-'^"'^''">l^<u\'^ (3) 
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with the generahsed index n — (/, j, k). The basis functions can 
have different widths (the k index) and different orders (the / and 
j indices). To deal with a spatially varying background a polyno- 
mial background term is present in Eq. [T]and is also included in 
the solution which is found by solving the resulting linear system 
(for details see A98 and AOO). The reference image is convolved 
with the kernel and subtracted from the search image. 

The default parameters for image subtraction used for the 
SVISS images are listed in Table [T] The basis functions de- 
scribed above used in ISIS 2.2 are three Gau ssian functions of 
orders 6, 4 and 3. Both llsrael etal] (|2007|) and lGossl & Riffesed 
(l2002i) find that the use of 3 Gaussian functions of these orders 
works well, we use these basis functions in our subtractions. We 
have, however, done some tests with different widths of the func- 
tions which is reported in section 15.31 The order of the poly- 
nomial func tions used for t he sp atial variation of the kernel in 
AOO is 2. IGossl & Riffesed (12002 ) use a zero order spatial vari- 
ability, but compute local kernels for subregions in the image. 
For the SVISS data set we subdivide the large 2k x 2.5k pix- 
els VIMOS frames into 9 subregions and use a first order spatial 
variation of the kernel. We have chosen to use kernel and stamp 
sizes scaled for the poorer seeing image (i.e. the search image), 
however we have not found these sizes to affect the subtraction 
quality substantially. A very small kernel size (smaller than the 
FWHM of the image) will likely cause the subtraction to fail. A 
too large kernel size will on the other hand include more noise 
in the fit and probably make the subtraction quality worse. In 
general, varying the parameters listed in Table [T] have consid- 
erably smaller effect on the resulting subtraction quality than the 
selection of stamps, data quality and host galaxy brightness (see 
section |43T l. 

In ISIS a number of sources present in both reference and 
search image are used as the data that goes into the kernel solv- 
ing procedure. A number of high S/N-sources are selected from 
the image and a small subimage (so called stamp) is obtained 
for each of them. The selected sources do not have to be point 
sources, but they have to be non-saturated and well below the 
non-linear response regime of the detector The stamps are se- 
lected by dividing the image into a number of rectangular sec- 
tions (as many as the wanted number of stamps). In each of these 
sections the brightest non-saturated pixel is found and the stamp 
is centred on that pixel. To check that no variable objects are 
selected, the fluxes inside the stamps in the two images are com- 
pared and the variable stamps are rejected. 

Our initial test runs showed some problems with the auto- 
matic stamp selection. The stamp selection is very sensitive to 
how saturated/bad pixels are handled in the data, if not treated 
correctly saturated objects could be selected as stamps. When 
simulating bright SNe we find that the simulated sources them- 
selves are often selected as stamps even though they are only 
present in one of the frames, and cause the subtraction qual- 
ity to deteriorate. This indicates that the rejection routine de- 
scribed above is not working for our data set with simulated SNe. 
However, often the number of SNe per image is expected to be 
quite low, thus the effect on real survey data is expected to be 
marginal. In our simulations, where a large number of artificial 
SNe were placed in the images the effect is clearly noticeable. 

To gain a better control over the stamp selection process 
we have modified ISIS to accept an external stamp list. This 
version of ISIS has previously been used for SN d etection in 
iMattila & Meikl^ (l2()01h : iMattila et all (|2004 12007|) with man- 
ual stamp selection. In this study the stamps are selected by run- 
ning SExtractor 2.4.4 (iBertin & Arnoutsll 19961) on the reference 
image, which allows us to reject stamps that are too faint and 



Table 1. Default subtraction parameters 



Parameter 


SVISS data set value 


Subregions 


9 


Total number of stamps {n^,) 


-450 


Size of the kernel (xFWHM) 


4.0 


Full size of the stamp (xFWHM) 


6.0 


Fitting order for the background 


2 


Degree of kernel spatial variation 


1 


Orders of the Gaussian basis functions 


6/4/3 



to place constraints on the stellarity of the sources. In this way 
we can also reject stamps that are suspected to contain a variable 
source and then rerun the subtractions. The stamp region list can 
also be compared with the bad pixel masks for both the refer- 
ence and the detection image to make sure that no stamp regions 
contain saturated pixels. 

In supernova surveys the sources selected as stamps are not 
necessarily point sources because the number of available bright 
point sources is simply too small. The majority of the auto- 
matically selected stamps will be extended sources, so if point 
sources are preferred the external stamp list method has to be 
used. In section 15.41 we have examined the effects of different 
selection criteria for the stamps and how the number of stamps 
affects the subtractions. 

2.3. Source detection 

Source detection is done in the subtracted frames using 
SExtractor (SE) to get an initial source list. The frames are 
also manually inspected to check for possible image defects 
and obvious spurious detections (e.g., sources close to satu- 
rated objects). The detection thresholds in SE are set by giving 
the minimum area of detected objects in pixels (MIN_AREA) 
and the minimum flux for a pixel to be considered detected 
(DET.THRESH). We have used MIN_AREA parameter corre- 
sponding to a circular area with diameter of roughly 1 xFWHM. 
The DET_THRES parameter is a multiplier to the measured 
pixel to pixel noise (rms) in the frame, this parameter has been 
adjusted to detect faint sources while still keeping the number of 
spurious detections down (see section 14.4.41 for a discussion on 
this). Convolving the images affects the noise and therefore also 
the optimal detection threshold to use, and since we are using 
other methods to ensure that spurious detections are excluded 
we have chosen to use a low detection threshold (1.5x the rms 
measured in the subtracted frame). 

In SE weight-maps in the form of variance maps can be used 
to correct the detections for varying noise levels over the image 
(due to e.g. vignetting and exposure time differences). This is es- 
pecially helpful in the edges of images. Using the SE option of 
WEIGHT.TYPE set to MAP .WEIGHT our exposure time maps 
are scaled to inverted variance maps within SE. The inclusion of 
a weight-map in the detection decreases the amount of spurious 
detections close to the edges of the subtracted images consider- 
ably. For each epoch weight-maps in the form of exposure time 
maps can be constructed (see|3]for a description on how this was 
done for the SVISS data set). For counts-per-second CCD im- 
ages where the instrumental noise is negligible (i.e. the noise is 
dominated by photon noise from the source or sky) the standard 
deviation of the individual pixel flux scales as the square root 
of the exposure time. The pixel noise in the subtracted image 
is approximately given by quadrature addition of the reference 
and search image noise (AL98). Convolving the reference im- 
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Fig. 1. Two example supernova candidates detected in the SVISS. Upper and lower panels show the two candidates with the leftmost 
panels showing the host galaxy in the reference image, the middle panels the host galaxy in the detection epoch and the rightmost 
panels the supernova candidate in the subtracted frame. The crosses in the leftmost panels mark the positions of the candidates. The 
candidate in the upper panel has an / band magnitude of 22.7 and in the lower panel 24.8. See text for more details. 



age slightly changes the noise characteristics of the subtracted 
image (discussed in AOO), but for constructing the combined 
weight-maps we have found it sufficient to use the unconvolved 
reference exposure time maps. The combined weight-maps are 
then computed as follows: 

= ^1^, (4) 

where Wref and Ws are the exposure time maps (weight-maps) 
for the reference image and the search image respectively and Wc 
is the weight-map for the subtracted image. Note that the weight 
maps are only used in the detection step in SE, photometry is 
done on the unweighted frame. 

All of the detections with valid photometry (in this case er- 
rors lower than 1.0 mag) are kept. The coordinates of the de- 
tected sources are all put into a master candidate list, this is then 
used for the photometry (described in the next section). The mas- 
ter candidate list can also be used to do multi-epoch and multi- 
band detection of the sources found in the search epoch image 
under consideration. The master candidate list will contain a 
large number of spurious detections (due to the low detection 
threshold used and the presence of subtraction residuals). These 
objects can be excluded by combining three different criteria; (i) 
too faint, the sources are fainter than a predefined limit (e.g. sig- 
nal to noise limit from simulations), (ii) only present in one filter, 
(iii) light curve not consistent with SN templates. 

2.4. Photometry of variable sources 

Photometry on the detected sources was done using the IRAF 
daophot package. PSF photometry was done on all the detected 
candidates using the task allstar. The PSF model used for 
each subtracted frame was constructed from the original worse 



seeing image (the search image) using the task ps£. In this task 
point sources are fitted with model profiles and the profile yield- 
ing the smallest residuals is chosen as the best PSF fit. The or- 
der of variability of the PSF model was set to zero, i.e. the PSF 
was not allowed to vary over the frame. Note that the model 
was constructed from the real data, thus the shape of the PSF 
is not symmetrical and does contain noise. As discussed earlier, 
often the PSF vary over the frames. However, using a higher or- 
der variability in the psf photometry did not result in any gain 
in the photometric accuracy. We also used aperture photometry 
with aperture corrections computed from the original worse see- 
ing image, using the IRAF task phot. This does give fairly good 
results, but seems to be more susceptible to residual flux from 
the background galaxies than the PSF photometry, thus giving 
larger errors (both statistical and systematic) in general. 

Using the pixel to pixel noise in the subtracted image will not 
give a correct estimate of the photometric errors, the pixels will 
be positively correlated due to the convolution of the reference 
image. Another problem in using the pixel to pixel noise esti- 
mate is that it will in general be computed in an annulus outside 
the aperture (or outside the size of the PSF model in the case 
of PSF photometry). The annulus region for most SNe will be 
outside the host galaxy, at least in the SVISS case where most 
host galaxies are just barely resolved. The noise estimate in this 
region can in principle be quite different from the noise at the 
SN position where residuals from the subtracted host galaxy will 
also be present. Both these effects result in the noise estimates 
given by daophot and phot underestimating the true photomet- 
ric errors severely. Instead we use the simulations described in 
section|4]to obtain reliable error estimates. By simulating SNe at 
different brightness in each epoch and finding the scatter in their 
measured magnitudes, we obtain an estimate of the true photo- 
metric error for each point on the SN light curve. 
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The limiting magnitude for a given epoch depends on the 
depths of both the reference and the search images. Since the 
noise in the subtracted image will be approximately equal to 
quadrature added noise of the images, the subtracted image will 
in general have a higher level of noise by a factor V2 and there- 
fore be ~ 0.4 magnitudes less deep than the individual images. 
For the SVISS data, this is true for the pixel-to-pixel noise. The 
subtracted frames do have on the order of V2 higher pixel rms 
(but slightly lower due to the convolution). However, we have 
found that the limiting magnitudes computed for the individual 
SVISS epochs based on pixel-to-pixel noise (m,,,,, „„5 in Table|2]l 
are not representative of the true depth as measured by by our 
simulation technique. We have measured the photometric scat- 
ter of the simulated SNe in the unsubtracted epoch images and 
obtained more conservative depth estimates (msim,sm in Table|2]). 

The reason why the photometric scatter estimate results in 
larger errors than the simple pixel-to-pixel estimate is a com- 
bination of several factors. We believe that the main factor is 
that our deep images are coming close to the confusion limit. 
In photometry the subtraction of a local background becomes 
very complicated when the background is allowed to vary over 
sizes comparable to the aperture size (or fitting radius for the 
PSF magnitudes). The background sources just at the confusion 
limit inside the photometric aperture/fitting region and back- 
ground annulus will vary from one simulated source to the an- 
other and can thus cause the scatter from the measured magni- 
tudes to go up. IncoiTect background subtraction in the presence 
of large scale background variations also seem to be a contribut- 
ing factor, although of lesser importance than the one described 
above. It is worth noting that subtracting the images actually 
serves to remove these two factors making the photometric scat- 
ter go down. In reality this effect is countered by the increase in 
pixel-to-pixel noise caused by the subtraction, see the discussion 
in section l5?n 



3. The SVISS data sample 

In this section we describe the data set that was used to in- 
vestigate the efficiency of the subtractio n method. The data 
were obtained with the VIMOS instrument (iLeFevre et al.l20()3h 
mounted on the ESO Very Large Telescope (UT3) at several 
epochs. The reference epoch was obtained in autumn 2003 and 
the search epochs autumn 2004. The VIMOS instrument has four 
CCDs, each 2kx2.4k pixels with a pixel scale of 0.205"/pxl, 
covering a total area of 4x56 sq. arcmin. The SVISS observa- 
tions were obtaine d in two fields, cove ring parts of the Chandra 
Deep Field-S outh (Gia cconi etaTIIIOO l) and the ELAIS-Sl field 
dLa Franca et a l. 2004), and in 5 broad band filters {U, B, V, R 
and /). The supernova search filters are R and /, with roughly 
twice the exposure time in / compared to R. Observations in 
these filters have been divided into 1 reference epoch (hence- 
forth epoch 0) and 7 search epochs. The time difference between 
epoch and epoch 1 is about one year and the search epochs are 
separated by a month. 

The data used in this work is a subsample of the SVISS, 
the simulations have been carried out in the deepest search band 
(the / band) and only on one of the four CCDs. The R band has 
only been used to determine the spurious detection ratio (see sec- 
tion l4.4.4l i. Table|2]contains an overview of the / band data used 
in the tests. It should be noted that the total number of frames in 
the reference epoch is ~ 80, we have only used the best frames to 
construct the reference frame in order to optimise the resulting 
seeing. The large number of frames available for making a ref- 



Table 2. Data sample overview 



Data set/Epocfi 


Exp. time (s) 






Seeing (") 


SVISS/Ref 


11520 


26.36 


26.1 


0.61 


SVISS/A 


11520 


26.64 


26.3 


0.68 


SVISS/B 


9600 


26.48 


26.1 


0.72 


SVISS/C 


14400 


26.62 


26.3 


0.64 


SVISS/D 


7680 


26.61 


26.2 


0.65 



t The limiting magnitudes are 3cr limits and have been measured for 
an aperture of diameter 4x the FWHM. 

X The limiting magnitudes are 3cr limits and have been estimated from 
the photometric scatter of simulated sources in the individual 
unsubtracted epoch images. 



erence image is not by design, most of the frames were obtained 
in a previous failed observation period (the VIMOS instrument 
was taken off-line due to technical problems in fall 2003). In this 
paper we do not discuss situations where the reference epoch 
has worse seeing than the search epochs. In the full SVISS data 
set we have some epochs with better seeing than the reference 
epoch. However, we have not used those epochs in the tests re- 
ported in this paper More details on the data and data reduction, 
along with source catalogues for the fields will be given in a 
forthcoming paper (Mencia-Trinchant et al. in prep). 

The data were reduced using a data reduction pipeline for 
SVISS developed by our team. The data were bias subtracted, 
flatfielded and scaled to correct for extinction. The images were 
also normalised to counts per second. A bad pixel mask was pro- 
duced for each of the images, containing vignetted (< 5% in 
most cases), saturated and cosmic ray flagged pixels. 

The VIMOS /-band suffers from rather severe effects of 
fringing, and care has to be taken to successfully remove these. 
We constructed a fringe map for each observation night by me- 
dian combining the non-aligned science frames and rejecting 
bright pixels by a standard sigma-rejection routine. The fringe 
map was then subtracted from each of the flatfielded science 
frames. All of the science frames were also manually inspected, 
some frames were removed because of vignetting and/or poor 
seeing (see Table |2]for the total number of frames used in each 
epoch). 

The good quality frames for each epoch were registered to 
a common pixel coordinate system using a simple shift and 
rotate transform (yielding a typical rms of ~ 0.1 pxls/0.02"). 
Finally, the registered frames were median combined. The bad 
pixel maps for each individual image were also registered and 
summed, yielding an exposure time map for each epoch. Some 
pixels (< 0.2% in the reference epoch) had an exposure time of 
0, i.e. none of the individual frames contained useful data for 
that pixel. These pixels were flagged in the combined science 
frames. 

Each of the combined epoch images were photometrically 
calibrated using secondary standard stars in the field. The cal- 
ibrated photometry for these stars were obtained from stacked 
one-night images from nights where standard star observations 
were available. 

For the reference image, a host galaxy catalogue was con- 
structed using SE. A total of 2660 galaxies were detected, re- 
jecting all sources with photometric errors larger than 0.1 mag 
and sources close to masked regions. For each of the galaxies a 
flux profile was calculated for four radii containing respectively 
20%, 40%, 60% and 80% (hereafter rio, r^, reo and rgo) of the 
total flux as measured by SE's AUTO photometry. 

The search epochs were aligned to the reference epoch fol- 
lowing the procedure outlined in section IZT] To obtain the nec- 
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Fig. 2. An example image from the SVISS data set containing a number of simulated supernovae. The field of view in the image 
is 1.9 'x 1.1 '(the total field of view of one VIMOS CCD is ~ 7 'x 8 ')• The top left frame is the reference image, the top right 
frame is the epoch 2 image with supernovae added (marked with with white crosses). The bottom left frame is the subtracted image 
(from run AO), the remaining sources are the added SNe. Some residuals are visible in the subtracted frame where no SNe have 
been added to the image (bottom right frame). 



essary accuracy for the registration we needed to use a transfor- 
mation that allowed separate shifts, rotations and scales for the 
X and y coordinates. The reference sources for registration were 
bright, point-like objects in the field (approximately 100 sources 
were used). The resulting rms of the geometrical transformations 
were smaller than 0. 1 pixels for all the epochs. 

The subtraction pipeline has a number of input parameters, 
which need to be adapted to the data set. Table [1] contain the list 
of parameters and the values for the SVISS data set. Since the 
VIMOS images are quite large we found it preferable to divide 
the images into subregions before the actual subtraction, which 
makes it possible to use a low degree of spatial variation in the 
kernel. The total number of stamps given in Table [T]is for the en- 
tire image. If the number of stamps was too small the subtraction 
quality degraded considerably. Since the number of stamps need 
to be tuned to the particular data set, we do not report any tests 
on varying this parameter The parameters related to the size of 
the stamp regions have not been found to affect the subtraction 
quality significantly. However, no sub-FWHM sizes have been 
considered. The parameters given in Table[T]have been kept con- 
stant in all the reported tests. 

Figure [T] show two supernova candidates detected in the 
/ band using the subtraction parameters discussed above. The 
bright SN candidate (upper panel in the figure) has an approxi- 
mate / band magnitude of 22.7 at peak and is situated inside a 
bright , tentatively elliptical, galaxy. The faint candidate (lower 
panel) has an approximate / band magnitude of 24.8 at detec- 
tion, and the explosion date for this candidate is not very well 
constrained. The host galaxy for this SN is fainter and more rep- 
resentative of average galaxies in our images. These detections 



demonstrate that the use of our subtraction pipeline enables us 
to find faint supernovae against the bright background emission 
from the host galaxies. 



4. Simulations 

4. 1 . The supernova point spread function model 

The construction of the artificial PSF model and the addition 
of SNe to image frames has been done with the IRAF package 
daophot (see also section l24b . 

For the SVISS data set approximately 40 bright star-like ob- 
jects were selected by visual inspection. About half of these ob- 
jects were rejected in the fitting procedure. The total number of 
point sources used for constructing the PSF model for each im- 
age was thus approximately 20. The best fitting profile was found 
to be an elliptical Moffat profile with 13-2.5. The FWHM of the 
models for different epochs/data sets can be found in Table|2l 

Supernovae of different brightness were simulated in the 
search images. In all test runs the SNe were added using the 
task addstar. The input model for this task was the best-fitting 
PSF model and the stars were added in the search image, i.e. 
the image with worse seeing. A total of 17 magnitude bins were 
used, the brightest bin roughly corresponds to the lOOcr limiting 
magnitude and the faintest bin somewhat fainter than the unsub- 
tracted image limiting magnitude (see Table |2]i. The magnitude 
bins used are: 22.0, 22.5, 23.0, 23.5, 24.0, 24.25, 24.5, 24.75, 
25.0, 25.25, 25.5, 25.75, 26.0, 26.25, 26.5, 26.75 and 27.0. 
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Run 


Stamp selection 


Host brightness 


Kernel widths 


O-rs/o-,-! 




df 




AO 


High SNR 


/ < 25 


(1.0,2.0,4.0) 


0.026/0.034 


1.00 


0.13 


0.825 


ASl 


Auto 


/ < 25 


(1.0,2.0,4.0) 


0.024/0.031 


2.97 


0.22 


0.848 


AS2 


Point 


/ < 25 


(1.0,2.0,4.0) 


0.024/0.034 


2.87 


0.12 


0.891 


AHl 


High SNRt 


/ < 22 


(1.0,2.0,4.0) 


0.027/0.077 


7.34 


0.47 


0.825 


AH2 


High SNR 


22 < / < 24 


(1.0,2.0,4.0) 


0.027/0.038 


1.76 


0.27 


0.825 


AWl 


High SNR 


/ < 25 


(0.2,0.4,0.8) 


0.027/0.033 


1.39 


0.13 


0.831 


AW2 


High SNR 


/ < 25 


(0.4,0.8,1.6) 


0.026/0.033 


0.924 


0.13 


0.821 


AW3 


High SNR 


/ < 25 


(0.6,1.2,2.4) 


0.026/0.032 


0.855 


0.10 


0.816 


AW4 


High SNR 


/ < 25 


(2.0,4.0,8.0) 


0.024/0.034 


1.88 


0.11 


0.834 


AW5 


High SNR 


/ < 25 


(4.0,8.0,16.0) 


0.027/0.035 


1.02 


0.14 


0.826 


BO 


High SNR 


/ < 25 


(1.0,2.0,4.0) 


0.023/0.035 


1.00 


-0.19 


0.667 


BSl 


Auto 


/ < 25 


(1.0,2.0,4.0) 


0.025/0.034 


1.25 


-0.32 


0.726 


BS2 


Point 


/ < 25 


(1.0,2.0,4.0) 


0.026/0.036 


1.34 


-0.29 


0.808 


BHl 


High SNRt 


/ < 22 


(1.0,2.0,4.0) 


0.024/0.084 


2.91 


-0.088 


0.667 


BH2 


High SNR 


22 < / < 24 


(1.0,2.0,4.0) 


0.024/0.041 


1.06 


-0.28 


0.667 


BWl 


High SNR 


/ < 25 


(0.2,0.4,0.8) 


0.024/0.034 


0.638 


-0.030 


0.648 


BW2 


High SNR 


/ < 25 


(0.4,0.8,1.6) 


0.023/0.033 


0.823 


-0.13 


0.659 


BW3 


High SNR 


/ < 25 


(0.6,1.2,2.4) 


0.023/0.033 


1.04 


-0.20 


0.672 


BW4 


High SNR 


/ < 25 


(2.0,4.0,8.0) 


0.024/0.035 


1.06 


-0.22 


0.677 


BW5 


High SNR 


/ < 25 


(4.0,8.0,16.0) 


0.026/0.046 


0.672 


-0.14 


0.661 


CO 


High SNR 


/ < 25 


(1.0,2.0,4.0) 


0.026/0.034 


1.00 


-0.21 


0.686 


CSl 


Auto 


/ < 25 


(1.0,2.0,4.0) 


0.027/0.035 


0.83 


-0.18 


0.674 


CS2 


Point 


/ < 25 


(1.0,2.0,4.0) 


0.032/0.041 


1.32 


-0.25 


0.808 


CHI 


High SNRt 


/ < 22 


(1.0,2.0,4.0) 


0.028/0.070 


2.92 


-0.69 


0.686 


CH2 


High SNR 


22 < / < 24 


(1.0,2.0,4.0) 


0.028/0.039 


1.40 


-0.51 


0.686 


CWl 


High SNR 


/ < 25 


(0.2,0.4,0.8) 


0.026/0.032 


0.921 


-0.20 


0.674 


CW2 


High SNR 


/ < 25 


(0.4,0.8,1.6) 


0.026/0.034 


1.19 


-0.25 


0.693 


CW3 


High SNR 


/ < 25 


(0.6,1.2,2.4) 


0.026/0.035 


0.949 


-0.20 


0.685 


CW4 


High SNR 


/ < 25 


(2.0,4.0,8.0) 


0.026/0.034 


1.06 


-0.20 


0.688 


CW5 


High SNR 


/ < 25 


(4.0,8.0,16.0) 


0.026/0.036 


1.18 


-0.21 


0.701 


DO 


High SNR 


/ < 25 


(1.0,2.0,4.0) 


0.025/0.033 


1.00 


-0.0053 


0.734 


DSl 


Auto 


/ < 25 


(1.0,2.0,4.0) 


0.025/0.034 


1.19 


0.048 


0.737 


DS2 


Point 


/ < 25 


(1.0,2.0,4.0) 


0.029/0.039 


1.19 


0.022 


0.832 


DHl 


High SNRt 


/ < 22 


(1.0,2.0,4.0) 


0.025/0.069 


4.00 


0.21 


0.734 


DH2 


High SNR 


22 < / < 24 


(1.0,2.0,4.0) 


0.027/0.038 


1.10 


0.099 


0.734 


DWl 


High SNR 


/ < 25 


(0.2,0.4,0.8) 


0.072/0.061 


2.15 


0.10 


0.821 


DW2 


High SNR 


/ < 25 


(0.4,0.8,1.6) 


0.030/0.036 


1.06 


0.033 


0.738 


DW3 


High SNR 


/ < 25 


(0.6,1.2,2.4) 


0.028/0.036 


1.18 


0.071 


0.738 


DW4 


High SNR 


/ < 25 


(2.0,4.0,8.0) 


0.029/0.040 


1.27 


0.079 


0.741 


DW5 


High SNR 


/ < 25 


(4.0,8.0,16.0) 


0.028/0.040 


1.66 


0.044 


0.749 



and CTri is the pixel to pixel noise in the isolated and inside galaxies positions, respectively, xiorm the normalised statistic for the run, df 
is the flux difference statistic and r,p is the spurious detection ratio (see section l4.4.1l for details), 
t Bright sources with simulated SNe have been manually removed from the stamp list. 



4.2. Placing the artificial SNe 

When adding artificial supernovae to the images care must be 
taken that the subtraction process itself is not affected. This puts 
some constraints on how many SNe can be placed in each im- 
age, and also where they are placed. Since ISIS relies on the 
use of stamps, it is important that there is a sufficient number of 
sources left without artificial supernovae. The image subtraction 
process should, however, be able to deal with the presence of a 
supernova inside a previously selected stamp (and reject it from 
the stamp list). In a real survey it is very unlikely that more than 
one supernova will be found in the same galaxy during the same 
survey period. The number of SNe per galaxy has therefore been 
limited to one. 

With these considerations in mind we are selecting host 
galaxies at random from a galaxy list, which, depending on the 
specific test run, can contain different samples of the full galaxy 
catalogue. The positions of the artificial sources inside the host 
galaxies are selected based on the flux profile for the specific 



host galaxy (in the general case, some test runs use a slightly 
different setup). Four position bins were used, where the in- and 
outside radii of the bins are given by: to r2o; /"2() to r4o; ^40 to 
W, '"60 to rgo. The angle with respect to the major axis of the 
galaxy and the radius within each bin were randomly selected. 
A fifth "position bin" was also considered, namely isolated po- 
sitions, i.e. positions in the image without any detected galaxy 
within a given isolation radius. These positions were found by 
masking the detected sources in the reference image, using a 5cr 
cutting threshold, dividing the resulting image into rectangular 
boxes with half-sides corresponding to the isolation radius. For 
the SVISS data an isolation radius of 20 pixels ~ 4"was used. 
The wanted number of supernova positions were then selected 
at random from this list of boxes and sources placed in the mid- 
dle of the boxes (one source per box). 

In Fig.|2]an example of an image with artificial SNe is shown, 
together with the subtracted frame. The lower right panel of 
the figure shows the subtracted frame without any simulated 



8 



J. Melinder et al.: Detection efficiency and photometry in supernova surveys 



sources. Typically, each subtracted frame will contain hundreds 
of residuals similar to the ones in this figure. 

4.3. Tests using the SVISS data set 

We have carried out a number of test runs with simulated super- 
novae added to the real images to study the effects on the sub- 
traction quality and SN detection. A summary of the properties 
of the different runs is presented in Table [3] for the SVISS data 
set. The tests have been done using four different epochs (the 
A, B, C and D designations, see Table |2] for basic properties of 
the different epochs). The reference runs for each epoch, AO-DO, 
are tests using an identical set of parameters. For each epoch we 
have then changed one parameter at a time to study its effect on 
the subtraction quality. Some parameters have been studied in 
many epochs, while some have only been studied in the epoch 2 
data. The same reference epoch (epoch 0) has been used for all 
the runs. 

The stamp selection has been done in three different ways 
(e.g. AO, ASl and AS2). In the "High SNR" option only the 
brightest sources in the field have been selected as stamps and in 
the "Auto" option the ISIS internal automatic stamp selection has 
been used. For the "Point" option, the most point-like sources 
have been selected as stamps. The total number of stamps in 
each of the three settings was ~ 450 over the entire frame (50 
per subregion). 

The tests on host galaxy brightness (i.e. 0, H1-H2) were done 
by selecting hosts from a magnitude limited galaxy list. In the 
H1-H2 runs the number of simulated SNe was a bit lower than 
in the reference runs, 100 instead of 266. For the "/ < 22" option 
the 100 brightest sources were selected as hosts, the mean bright- 
ness of these galaxies was / - 20.0. For the "22 < / < 24" option 
100 sources were selected at random from a magnitude limited 
catalogue, the mean brightness of these sources was / - 23.2. 
The mean brightness of galaxies in the full host galaxy cata- 
logue was / - 23.7. In the case of the brightest host galaxies, 
some stamps had to be removed manually as they were not re- 
jected by ISIS despite the addition of a bright SN to the stamp 
region. Having SN host galaxies as stamps affected the subtrac- 
tion quality strongly for the brightest supernova bins. 

In the W1-W5 runs, the width of the Gaussian basis func- 
tions that constitute the fitting kernel have been changed. The 
default ISIS values for the widths (expressed as the Gaussian cr 
parameter) are (1.0,2.0,4.0). In the test runs we are using both 
narrower and broader widths. The optimal kernel should contain 
a basis function with width that when convolved with the refer- 
ence image will produce an image with the same seeing as the 
search image. For the SVISS images (see Table |2]i the optimal 
kernel should thus contain a basis function with cr ~ 0.4 - 0.8. 

4.4. Test diagnostics 

4.4.1. Pixel statistics in subtracted images 

In the AL98 and AOO papers pixel statistics of the subtracted 
images are used as a measure of the subtraction quality. The the- 
oretical pixel noise distribution (hereafter PND) is compared to 
the measured distribution in the subtracted frames. The images 
used in these subtractions are simulated images with a purely 
Gaussian PND. The theoretical distribution is thus quite sim- 
ple to calculate, even though convolution of the reference im- 
age makes it somewhat more complicated. In our tests we have 
used real data images with many sources, both point sources 
and extended objects, the PNDs for these images are harder to 




m . 



Fig. 4. Detection efficiency for isolated positions and different 
galactocentric distances inside the host galaxies. The plot was 
constructed from the AO run. (This figure is available in colour in 
the electronic version of the article) 

model. We have opted to take a more observationally orientated 
approach. Instead of comparing the PND in the subtracted im- 
age with a theoretical distribution, we compare the PND of iso- 
lated positions in the subtracted image to that of positions inside 
galaxies (see Fig.[3]for an example of this). The PND for isolated 
positions is thus considered to be the optimal one, the closer the 
PND of the inside galaxy positions is to the isolated case, the 
more optimal the subtraction. 

The pixel lists used in the statistics contain all pixels in- 
side a 3" radius from all of the simulated SN positions, i.e. the 
pixel lists are obtained from a list of apertures, measured in an 
"empty" (without simulated SNe) subtracted frame. The posi- 
tions for the isolated case are taken from the fifth (isolated) po- 
sition bin and the inside galaxies positions from the rl position 
bin. The same number of apertures have been used for both pixel 
lists, so that the total number of pixels is the same within inter- 
polation errors. 

For each of the test runs we have calculated three different 
statistics based on comparing the PNDs: (i) the pixel to pixel 
noise inside the apertures of the inside positions and the isolated 
positions, ctii and cTi-s', (ii) the reduced;i'^ statistic computed 
from comparing the PNDs; (iii) the flux offset statistic, df. 

The flux offset statistic is obtained by numerical integration 
of the PND difference over the counts per pixel and dividing by 
the number of apertures used. The flux offset is defined as: 

df^—-{\ -{Pn{x)-Piso(x))dx 
riaiAJii 

+ j^^(Pri(x)-PM)dx'^, (5) 

where nap is the number of apertures, // and ul are the upper 
and lower limit of the counts per pixel, respectively, x counts 
per pixel and P^i /iso the number of pixels in the interval dx for 
a given value of x. The rl and iso designations refer to the two 
position bins considered (see above). This statistic is in princi- 
ple the mean counts per aperture with the background, defined 
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Fig. 3. Pixel noise distributions (PNDs) for two of the search epochs (upper panel shows search epoch A, lower epoch B). Leftmost 
panels show the reference test runs (using high-SNR stamps and completely random host galaxies), middle panels show the S 1 runs 
(using ISIS automatic stamp selection) and rightmost panels show HI runs (bright host galaxies). All distributions, except the bright 
host galaxy test, contain the same total number of pixels (for details see Table[T]and section l4.4.1b . The dashed (red) line show the 
difference of the inside and outside pixel distributions for each case. (This figure is available in colour in the electronic version of 
the article) 



as the counts inside the isolated position apertures, removed. The 
numerical integration of the PND difference enables us to con- 
trol the lower and higher bounds as well as making sure that the 
correct background, i.e. the counts within isolated position aper- 
tures, is removed. 

4.4.2. Detection efficiency 

For each of the magnitude bins a number of detected sources is 
obtained from the SE runs. A detected source is defined as hav- 
ing a signal-to-noise ratio (SNR) of more than one, "sources" 
detected by SE can have very low significance and include many 
spurious detections. The input simulation coordinate list is then 
used as an ASSOC list in SE (i.e. matched in image coordinates) 
to make sure that only simulated objects are counted in the de- 
tection. In this way we can study the detection efficiency for a 
certain set of parameters in real data without including spuri- 
ous detections from residual sources. This estimate of the sur- 
vey depth is useful when other methods can be used to dispose 
spurious detections, e.g. by using multiple bands or epochs. For 
single epoch measurements greater care has to be taken not to 
include spurious objects. 

The detection efficiency is then calculated for each magni- 
tude bin as the number of detected sources divided by the to- 
tal number of simulated objects in the magnitude bin. However, 
these will not be the true detection efficiencies since the SNR es- 
timated by SE is not accurate due to the correlated pixel to pixel 
noise in the subtracted images. Such detection efficiencies would 
thus be unreliable to use in the SN rate calculations, and in this 



paper we use them only for comparison between the different 
subtraction parameters and epochs. Due to the way we handle 
detections in the images with simulated SNe (only accepting de- 
tections at positions where SNe were placed), the signal to noise 
cut used does not affect the results discussed in this paper For 
the SN rate calculations we will use detection efficiencies based 
on signal to noise cuts found from doing PSF photometry of sim- 
ulated SNe (using the photometric scatter method described be- 
low). 



4.4.3. Pliotometric accuracy 

The photometry of the supernova candidates is done in the sub- 
tracted images using PSF photometry (see section |24l ). The ac- 
curacy of this method has been assessed by subtracting the input 
simulation magnitudes from the measured magnitudes for each 
of the simulated sources and then calculating the mean, median 
and standard deviation of the magnitude differences. The magni- 
tudes have also been converted to flux (in counts/sec) using the 
magnitude zeropoints. The difference in flux and magnitude ver- 
sus input simulation magnitude have then been plotted for the 
different cases. It is important to note that each of the magnitude 
bins is independent of each other, the host galaxy and simulated 
SN position are selected at random for each of the bins. Any 
offset in the the accuracy plots thus means that the entire im- 
age/method suffers from a systematic error, rather than effects 
from selecting a "bad" sample. Tests have also been done us- 
ing aperture photometry with aperture correction, this gives very 
similar results but with a larger scatter. 
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Fig. 5. Photometric accuracy for SNe at different galactocentric distances as a function of simulation magnitude (msim). The plot 
was constructed from the AO run. The upper left panel shows the median pixel flux difference (in counts/sec) between the input 
simulation flux (Fsim) and the measured PSF flux (Fpsp), and the upper right panel shows the median magnitude difference. Error 
bars in the flux difference plot show the standard error of mean of the difference. Note that the detections become incomplete at 
magnitudes > 26, this causes the apparent failure of photometry at those faint magnitudes. The lower left panel shows the signal 
to noise versus the simulation magnitude and the lower right panel the measured scatter in each simulation bin. The signal to noise 
ratios have been estimated from the measured scatter (SNR ~ 1.0857 x -^f—)- The dashed horizontal line marks the SNR - 3 limit. 

(This figure is available in colour in the electronic version of the article) 



The standard deviation of detected SNe in each magnitude 
bin is used as an estimate of the photometric scatter As dis- 
cussed in section |24l the error estimates from daophot and phot 
will not give accurate photometric errors. The standard deviation 
in each simulated magnitude bin should provide a more accurate 
estimate for the statistical photometric errors. The SNR for a 
given magnitude is computed from the magnitude scatter using 
the formula SNR - 1.0857 x 

A problem with estimating photometric accuracy in this way 
is that the detected sample will be incomplete at the faintest 
magnitudes. Only the very brightest sources in a given magni- 
tude bin (i.e. sources that happen to be located at the position 
of a noise peak or residual flux peak) will be detected at the 
faint limit. This can be seen in the photometric accuracy plots, 
e.g. Fig. |5] where the flux difference becomes positive (indicat- 
ing that the flux is overestimated) in the faintest magnitude bins. 
Conclusions should only be drawn based on the complete bins. 

4.4.4. Spurious detections 

Whenever source detection codes are used to do automatic de- 
tections there is a risk that some of the detected sources are in 
fact not real, but rather noise peaks that are mistaken for real ob- 
jects. The number of spurious detections is affected by the detec- 



tion parameters, a very conservative setting of these parameters 
can result in few spurious detections. However, the problem with 
using conservative detection parameters is that real, but faint, 
sources might be missed. In galaxy surveys it is important that 
the number of spurious detections is known and, in some cases, 
controlled. Different methods can be used to do this, one is to 
run the detection again with the same parameters on the inverted 
image. If the noise is Gaussian the amount of spurious detections 
is expected to be similar in both the images. 

For the subtracted images used in this work this method is 
not suitable. Many of the spurious detections will not be re- 
lated to the Gaussian noise in the isolated parts of the image, 
but rather associated with residuals from galaxies that were not 
fully subtracted. The "noise" in these residuals is not necessarily 
Gaussian nor symmetric, thus the inverted image is not a good 
tool to use. 

We instead define a spurious detection in the subtracted im- 
ages as a source that is detected in the test search epoch / band 
images, but not detected in the R band images for the same 
epoch. The spurious detection ratio is thus defined as: 



where «/ is the number of detected sources in the / band image, 
which is the filter used for all of the tests. njcR is the number 
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Fig. 7. Photometric accuracy for the different search epochs. Panels are the same as in Fig.|5] The legend shows the symbols used 
for the epochs. (This figure is available in colour in the electronic version of the article) 
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Fig. 6. Detection efficiencies for isolated and inside galaxies po- 
sitions at the four different search epochs. (This figure is available 
in colour in the electronic version of the article) 



5. Results 

Our tests show that the cr,.i/cr,.5, x^^^j and df statistics can be 
used as measures on how successful the difference imaging was 
in subtracting galaxies from positions where artificial SNe were 
added. We have also normalised the xj.^^ statistic for each test 
run by the value found for the reference run (AO, BO, CO and 
DO). The normalised x^^^ makes it easier to compare the test 
runs within each epoch, but can not be used to compare epochs 
to each other Table [3] contains the resulting statistics for the test 
runs, and Fig. [3] shows the pixel noise distributions for some of 
them. 

Some general remarks can be made on the relationship be- 
tween the different statistical measures. By definition will be 
a measure of the total difference of the two distributions. There 
are two main components to this total difference, one is the dif- 
ference in widths of the two distributions (characterised by the 
standard deviations, (Xri/crrs). The df statistic is mainly sensitive 
to possible systematic flux offsets, which is the other main com- 
ponent in the ^^^.i^, and gives an indication on whether the offset 
is positive or negative. In Table[3]we see that, as expected, a large 
CTi-i and/or a large df (negative or positive) results in a large 
This is also evident from looking at the pixel distribution plots 
in Fig. [3] 



of sources detected in both / and R. Source detection using SE 
is done in "empty" subtracted images for the different param- 
eter sets. For this test the detected sources are from the entire 
image, not only from the simulated positions considered for the 
detection efficiency tests. 



5. 1 . Detection and photometry of supernovae with faint host 
galaxies 

To test the effects of the supernova galactocentric distance we 
placed the simulated sources as described in section l4~2l The 
resulting detection efficiency and photometric accuracy for the 
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different position bins are shown in Figs. |4] and |5] While the po- 
sition tests have been done in all the test epochs, these figures 
show only run AO. The four position bins inside each galaxy 
were merged into two, one region containing 40 % of the flux 
and the other one containing 40 to 80 % of the flux. We find no 
difference within the eiTors between the inside position bins, nei- 
ther for detection efficiency nor photometric accuracy. The small 
differences visible in the figures is smaller than the scatter for 
each magnitude bin. The scatter in the detection efficiency plot 
is hard to calculate, but the fact that for values for the two posi- 
tions switch order between magnitude bins can be interpreted as 
a scatter effect. This result holds true for all of the test runs. 

The resulting efficiency and photometric accuracy for the 
isolated positions differs significantly from the ones for inside 
positions. Assuming the isolated SN positions to be an ideal case 
and comparing the photometric accuracy for these positions with 
the inside ones, we find that there is a systematic error in the pho- 
tometry for the inside positions. Further testing showed that this 
systematic error is in fact a flux offset, the difference between 
simulated and measured flux is the same in all of the magnitude 
bins (shown in the left panels of Fig. |5] and IT). In the case of 
run AO the flux offset is positive, positions inside galaxies have 
too much flux remaining after subtraction. This is also visible in 
the detection efficiency plot, the efficiency for inside positions is 
erroneously boosted by the extra flux, giving an efficiency bet- 
ter than the isolated case, which is clearly incorrect. The eiTor 
bars given in the flux difference panels show the standard error 
of mean in each bin. The difference is nonzero (within 3-5 cr) 
and is consistent with a constant flux offset in all of the magni- 
tude bins. Using the mean flux of the faint host galaxy sample 
and the measured flux offset we can calculate the mean fraction 
of host flux that remains in the subtracted frame. For the A, B, 
C and D epochs the fraction is 1.9 %, 3.0 %, 2.5 % and 5.3 %, 
respectively. 

The subtraction has succeeded to varying degree in the four 
different search epochs, the detection efficiency and photomet- 
ric accuracy for the reference runs of each epoch is shown in 
Figs. |6] and |7] as well as in Table |3] The x^^j statistic for the 
four reference runs are in fact very different, 1 1 .9, 46.2, 32.7 and 
17.3, respectively, for AO-DO. This is not surprising since the 
data from the different epochs are of somewhat differing qual- 
ity (see Table |2]i. The flux offset has been detected in all of the 
test epochs to varying degree, in epoch A it is a positive offset 
and in B, C and D a negative one. This is clearly visible both 
in the trend for the df statistic and in the photometric accuracy 
plot (Fig.|6]l. In all of the epochs, a flux offset also results in an 
artificial boost/suppression of the detection efficiency. 

The lower panels of Fig. |5] and [T] show the SNR and stan- 
dard deviation of the measured magnitudes in each bin for the 
different positions and epochs. Contrary to the standard error of 
mean, the scatter is an estimate of the photometric uncertainty 
for individual supernovae. As expected, the scatter is consider- 
ably smaller in the isolated positions compared to the inside po- 
sitions. The scatter and SNR for the three epochs are very sim- 
ilar, but the detection efficiency is different. This shows that the 
detection efficiency is affected by the systematic flux offset. The 
SNR for a given simulated magnitude bin is higher for the iso- 
lated positions than for the inside positions (lower left panel of 
Fig-El. This is of course as expected since the presence of a host 
galaxy will affect the photometry of the SN. The estimates of 
SNR and scatter are unreliable for the faintest magnitude bins 
since the sample becomes incomplete in these bins (e.g., see 
Fig.m. A useful estimate of the depth can be found by extrapo- 
lating from the trustworthy data points where the detections are 



100 % complete. The 3cr detection limit for the isolated position 
is approximately nij - 26.1. This is comparable to the 3cr lim- 
iting magnitude for the individual frames (see Table |2]i. The in- 
crease in noise from subtracting the images which should cause 
the limiting magnitudes in the subtracted images to be on the or- 
der of 0.4 magnitudes lower is not seen. This is due to the effects 
discussed in section |Z4l where the subtraction actually serves to 
remove the problematic background, this effect is roughly sim- 
ilar (by coincidence) to the increase in pixel to pixel noise. The 
SNR versus magnitude plots for the four different epoch (lower 
left panel of Fig. |7]i show that the B epoch is considerably less 
deep than the other three epochs. This is in agreement with the 
limiting magnitudes given in Table |2l where epoch B is indeed 
the shallowest of the four epochs. The 50 % completeness limit 
roughly corresponds to the 3cr limiting magnitude obtained from 
the SNR plots, but it should be noted that the flux offsets have a 
fairly strong effect on the detection efficiencies. The number of 
spurious detections, characterised by r^p, also reflects this, the 
number is significantly higher in the A epoch (which has a posi- 
tive flux offset) than in the epochs with negative offsets. 

5.2. Detection and photometry of supernovae with bright 
host galaxies 

Properties of the runs with different host galaxy characteristics 
can be found in Table[3](Hl and H2 runs for the 4 different search 
epochs). Figs. [8] and |9] show the detection efficiency and photo- 
metric accuracy in the reference, HI and H2 runs for one of the 
epochs (epoch A). 

The detection efficiency plot clearly shows that the efficiency 
becomes worse when the host galaxy is considerably brighter 
than the simulated supernovae. This could possibly be explained 
by the Poisson noise from the host galaxy eventually becoming 
dominating when fainter SNe are detected. Assuming that the 
noise for a source consists only of Poisson noise from the host 
galaxy, i.e. SNR ~ /^sn/ V^kKt, we find that the SN magnitude 
for a 90 % detection efficiency coiTesponds to a signal-to-noise 
ratio of ^ 100 for the bright host galaxy group. Thus it seems 
unlikely that the host galaxy Poisson noise affects the detection 
probability nor the photometry for the simulated SNe. Also, the 
drop at the detection limit is considerably slower for the bright 
and intermediate brightness hosts. These two effects can be ex- 
plained as subtraction residuals from the host galaxies being de- 
tected as SNe. Visual inspection of the subtracted frames also 
show that almost all of the selected bright host galaxies show 
significant residuals. The pixel-to-pixel noise in apertures placed 
on the inside galaxies position (cr,-i given in Table [3]) is also 
increased substantially when considering brighter hosts. This 
is due to the subtraction residuals in these positions since the 
Poisson noise from nearly all but the very brightest galaxies is 
negligible in our deep and sky noise dominated images. 

The photometric accuracy is also affected by the brightness 
of the host galaxies. The constant flux offset discussed in the 
previous section is no longer seen for the two populations of 
brighter host galaxies (see Fig. |9]l. In the very brightest galaxy 
population we see that the flux and magnitude differences are 
very large. There seems to be some kind of systematic offset 
inside the bright galaxies but it is no longer constant in the dif- 
ferent magnitude bins. For this host sample, the magnitude bins 
are not independent of each other, the host galaxies are the same 
in all the bins, i.e. the 100 brightest galaxies in the field. The flux 
difference is in this case not independent of the SN flux. For the 
very faintest bins the difference is positive (which also the df 
statistic for run AHl in Table [3] indicates), while for the very 
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Fig. 8. Detection efficiency for the runs with different host 
galaxy brightness for a single epoch (AO, AH1-AH2). The mean 
/ band magnitudes for the three host galaxy populations are 20.0, 
23.2 and 23.7. (This figure is available in colour in the electronic 
version of the article) 

brightest simulation bins the difference is negative. However, 
it should be noted that the scatter is very large for the bright 
hosts. The SNR plot for different host galaxy brightness (lower 
left panel of Fig.|9]l shows that the 3 cr limiting magnitudes for 
the bright and intermediate host galaxy samples are a lot worse 
than for the faint galaxy sample. 

A possible explanation for this effect could be that the pho- 
tometry fails in the presence of powerful negative residuals, giv- 
ing unpredictable results for the bright simulation bins. The in- 
termediate brightness hosts are somewhere in between the faint 
and bright hosts, for the epoch A case they do not show any offset 
at all for any of the magnitude bins. We cannot find any simple 
relation between the host brightness and the found photometric 
difference, the epochs are affected differently, c.f. AH-DH runs 
in Table [3] (Fig. [8] shows only epoch A). Since the subtraction 
residuals at the bright galaxy locations vary considerably from 
epoch to epoch, this is expected if the residuals is the source of 
the found photometric errors for bright host galaxies. 

It should be noted that the bright host runs were carried out 
with a slightly modified stamp list. We found that bright galaxies 
with simulated supernovae were not correctly rejected by ISIS 
in the subtraction. The subtraction quality was therefore greatly 
degraded for the brighter magnitude bins. We manually removed 
the bright host galaxies with simulated SNe from the stamp list. 
The stamp rejection is further discussed in section l54l 

5.3. Effects of convolution kernel properties 

Detailed testing on subtraction parameter s in OIS me t hods h ave 
already been reported by, e.g., AOO and llsrael et al.l (l2007h . In 
these papers the detailed testing is done on fully simulated data 
where noise and subtraction residuals can be controlled. We have 
done some tests on varying the subtraction parameters in our real 
data, i.e., observed images with simulated SNe. In general, the 
noise and residuals from the host galaxies dominate over any 
effects resulting from subtraction parameter changes. We have 



therefore chosen to use the default ISIS 2.2 parameters for all 
parameters but the convolution kernel width (see TablelTJ. In the 
reference runs (the runs) the default kernel widths have been 

used; 

llsrael et akl (l2007h find that the widths of the basis func- 
tions of the kernel affect the subtraction result quite strongly. 
We have tested this in our real data and the results are shown in 
Table |3] The tests show that the effect of changing the widths 
is small for most of the tested width combinations. The xlorm 
statistic indicate that the widest (W4-W5) kernel widths are the 
worst. The best kernel width parameter setting is different for 
the four epochs but in many cases kernel width combinations 
that are more narrow than the default give a better result. We 
have also calculated the detection efficiencies for several kernel 
width combinations. The and W1-W3 runs all have very simi- 
lar efficiencies, but for the widest width combinations the depth 
goes down considerably. This is consistent with the noise in the 
subtracted image (ctji) being ~ 4 times higher when using the 
widest basis functions (W4-W5) than in the reference case. In 
one case (BW5) the Xmnn statistic indicates that this is the best 
setting for that particular epoch, but the noise in the subtracted 
frame is high and the detection efficiency is very poor. For our 
data set, where the seeing differences are quite small, it is likely 
that the lack of at least one narrow basis function makes the con- 
volution kernel suboptimal. In principle one would expect the 
pixel to pixel noise to go down when using a wider convolution 
kernel, but it seems that the higher subtraction residuals in reality 
increase the noise substantially in isolated positions and slightly 
in the non-isolated positions. 

5.4. Stamp selection 

The different stamp selection methods are compared using the 
subtraction quality diagnostic statistics in Table [T] The tests 
show that the "High SNR" option works best, both in terms of 
overall subtraction quality, as estimated by the Xmmn ^^'^ di- 
agnostics, and in terms of spurious detections. The "Point" op- 
tion gives the worst results for all of the tested epochs. The 
"Auto" option gives differing results, which is not surprising 
since the stamp list will be remade for each epoch, in contrast 
to the two other methods where roughly the same stamps are 
used for all test epochs. In general, this method seems to give 
slightly worse results than the "High SNR" one. 

It is also worth noting that the stamp rejection routine does 
not seem to work properly. Galaxies with SNe (bright or faint) 
are not always rejected in the stamp selection by ISIS. This is 
a problem for the automatic stamp selection method where the 
subtraction quality can go down considerably. With the other 
methods the stamp list can always be manually edited before the 
subtraction and the non-suitable stamps removed from the list. 

5.5. Estimating photometric errors for SN light curves 

Running the detection pipeline on simulated supernovae can pro- 
vide useful estimates of possible systematic offsets as well as 
photometric errors. Figure [10] shows the light curve of a type 
la supernova at redshift 0.5 (light curve model constructed us- 
ing the templates provided by R NugenQ). The measurements 
were done in the subtracted images using PSF fitting photometry 
(as discussed in section l24b . Using the simulations reported in 
section ISTI we have corrected the light curve for the systematic 
flux offset and estimated the uncertainty for each epoch mea- 

' jhttp://supemova.lbLgov/~nugent/nugent_templates.html| 
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Fig. 9. Photometric accuracy for runs with different host galaxy brightness for a single epoch (AO, AH 1-AH2). Panels are the same 
as in Fig.|5] The legends show the symbols used for the different host galaxy tests. (This figure is available in colour in the electronic 
version of the article) 



surement. Since the pixel to pixel noise in the subtracted images 
is correlated the magnitude errors reported by daophot will not 
be reliable (they will underestimate the errors). The uncertainties 
are instead estimated as the scatter of the simulated magnitude 
bins (see figure|5] interpolating when necessary). The systematic 
offset and the scatter will be different for each epoch, thus sim- 
ulations need to be carried out in each epoch. These realistically 
estimated errors are on the order of 2 - 5 times higher than the 
standard output from daophot. 

The correction for the systematic flux offset is calculated as 
a mean quantity for a large number of SNe. For the case shown 
in Fig.[TO]the corrected magnitudes is a better match to the sim- 
ulated light curve than the uncorrected ones, the mean residual 
for the uncorrected and corrected curve is 0.175 and 0.0775, re- 
spectively. It is evident that the small flux offset is negligible for 
the brighter points on the light curve and it is mainly the faintest 
point that is affected by the correction (also the one giving rise 
to most of the change in mean residual). It should also be noted 
that once the faintest data point has been corrected it is consistent 
(i.e. within 1 cr) with the simulated light curve. 

6. Discussion 

6. 1 . Overall pipeline performance 

We have studied the detection efficiency and photometric ac- 
curacy of our SN detection pipeline using real data from the 
SVISS. However, these results can likely be applied to any SN 
ground based imaging survey using ISIS subtracted images for 
detection and/or photometry. Previous tests of image subtrac- 



tion software have been done mainly on artificial images (AL98, 
AOO, Gossl & Riffeser 2002; Israel et al. 2007) but also on stel- 
lar fields (AL98). To our knowledge the work presented in this 
paper is the first detailed study of subtraction quality in real 
supernova survey data. Detection efficiency studies in super- 
nova survey data are regularly carried out to be able to find the 
supe rnova rates (e.g. ^Cappellaro et al. 2005; Harris et al.l 120051 ; 
Mi knaitis et a n i2007h . but this testing is seldom reported in de- 
tail. 

We have found a systematic flux offset affecting the photom- 
etry in subtracted images. In general, with host galaxies chosen 
at random in the field, the offset is not dependent on the SN 
brightness, the offset is constant for different simulated SN mag- 
nitude bins. The offset is only found for SNe placed inside host 
galaxies whereas in isolated positions no offset is found. The off- 
sets change between the different epochs, in some cases being 
positive and in others negative. Pixel statistics for pixels inside 
apertures centred on the host galaxies also show that there is a 
shift of the median pixel value in the same direction as the flux 
difference plots indicate and as expected, a broadening of the 
distribution. We interpret this as a small flux through/plateau at 
the locations of galaxies, reaching as far out as to the 80% flux 
radius. Currently we do not have an explanation for this effect, 
it cannot be solved by changing the convolution kernel sizes or 
stamp sizes. Further testing is needed to pinpoint if the system- 
atic offset is caused by properties of the dataset or if it is in fact 
an internal effect of the image subtraction method. It should be 
noted that this is a small offset, the systematic error on the mag- 
nitudes of SNe is smaller for brighter sources (for the brightest 
SNe considered in this paper, my - 22.0, the offset is ~ 0.005) 
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Fig. 10. "Observed" light curve of a simulated type la supernova 
at z = 0.5. Solid line is the model light curve in the VIMOS 
I band from peak magnitude to ~ 120 days (in the observers 
frame). Open (black) squares show the measured SN magni- 
tudes in the subtracted frames and solid (blue) hexagons show 
the corrected magnitudes (see text for details). The error bars 
show the scatter measured by the detailed simulations reported 
in section ISTl and are shown as ±lcr error bars to the corrected 
magnitudes. For the brightest magnitudes the errors are smaller 
than or roughly the size of the corrected points and thus not visi- 
ble in the plot. (This figure is available in colour in the electronic 
version of the article) 



and in general significantly smaller than the photometric scatter 
for an individual source. However, in a survey of multiple ob- 
jects, and particularly fainter objects, this systematic effect will 
likely affect the result unless dealt with. 

A simple estimate of the photometric uncertainties using the 
pixel to pixel noise in the local background underestimates the 
errors when doing photometry on SNe in subtracted images. To 
get an accurate noise estimate we have instead computed the 
photometric scatter for each simulated magnitude bin. This scat- 
ter has also been used to construct signal to noise ratios, using 
this method the SNR is not affected by the flux offset and pro- 
vid es a simple way of rejecting less trustworthy observations. 

[Miknaitis et al. (2007) use simulations to check whether im- 
age registration and subtraction affect their SN photometry. They 
find no systematic offsets to a SNR of 5, but this requires an ac- 
curate SN position obtained from multiple points on the light 
curve. They find the photometric scatter to be on the order of 20 
% larger than what their noise maps indicate (but the noise maps 
already include the effects of re-binning and convolution). This 
seems considerably better than the accuracy we find using our 
pipeline and simulations on the SVISS data. The method used 
by these authors is quite different from the one used in this work, 
therefore it is hard to compare the exact figures. We estimate the 
SNRs differently by using the scatter found from simulations as 
our noise. It s hould also be noted tha t the SN positions in the im- 
ages used bv iMiknaitis et al.l ( l2007h are chosen at random (pri- 
vate correspondence G. Miknaitis), thus their results will not be 
directly comparable to our inside host galaxies simulations. The 
main difference is that we find that subtraction residuals from 



the host galaxies can have a noticeable effect on the photometric 
accuracy of the survey. This effect does not seem to be present 
in the work of Miknaitis et _aL, (2007) . 

The presence of spurious detections is a problem for all su- 
pernova surveys, and in particula r for surveys with sp ectroscopic 
follow-up (see the discussion in Bailey et al.l 120071) . For multi- 
band, multi-epoch imaging surveys (like the SVISS) the spurious 
detections is less of a problem, but it is still desired to keep the 
number down. The bright residuals found in less successful im- 
age subtractions can easily be mistaken for real varying sources 
when using automatic source detection. The number of spurious 
detections in our tests is large, partly because we use a low detec- 
tion threshold, but in general we find that a single epoch, single 
filter run results in a spurious detection rate of around 80% (i.e. 
1 out of 5 are real sources) when the detection threshold is set 
to find objects brighter than the 5cr limiting magnitude. For the 
SVISS dataset used in this paper the total number of detected 
sources for a given epoch is on the order of hundreds, making 
the number of spurious detections quite large. We find that the 
use of a controlled and bright stamp list makes the number of 
spurious detections go down by about 10%. It should be noted 
that residuals at the locations of bright galaxies will be problem- 
atic even in a multi-band survey since the residuals often appear 
at the same positions in all filters. 



6.2. Host galaxy brightness and supernova position 

We find that the typical brightness of the host galaxies does af- 
fect the detection efficiencies and photometric accuracy. This is 
of c ourse as expected and also confirmed by other SN surveys 
(e.g. Neill et al. 2006). However, both detection efficiency and 
photometric accuracy becomes gradually worse when consider- 
ing host galaxies of increasing brightness. The effect is mainly 
noticeable for the very brightest host galaxies (i.e. largest con- 
trast between the SN and the host). The fact that the photometry 
is affected by host brightness can in principle cause the system- 
atic offset found above to vary with the SN brightness, since the 
SN and host brightness are correlated through the host galaxy 
distance. In projects where the supernova rate is studied, the ef- 
fects of host brightness most likely have an impact on the detec- 
tion efficiencies at low redshifts, where the host galaxies can be 
very bright, as well as in bright galaxies at high redshifts where 
the supernovae will be very faint. Whether this effect will be no- 
ticeable depends on the characteristics of the galaxy population, 
and simulations of the type reported in this work can success- 
fully quantify the effect for galaxies of different brightness. This 
information can then be incorporated into the overall detection 
efficiencies discussed above. The SNe detections in the SVISS 
data are likely to be found in faint galaxies and thus in the gen- 
eral case not affected strongly by the host background flux. 

The effect of the SN galactocentric distance in the host 
galaxy is difficult to study in the SVISS data set, since most of 
the galaxies have small angular sizes, and in many cases are not 
even resolved. Our results indicate that the flux offset is slightly 
worse in the inner parts (40% in flux units) of the galaxies than 
in the outer (40-80%) regions. Therefore, we cannot draw any 
conclusions on how the detection efficiencies are affected by the 
galactocentric distances. However, if the detection efficiency is 
dependent on radial distance from the host centre, the radial dis- 
tribution function of SNe needs to be known to find a proper 
correction factor for the total efficiency. 
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6.3. Image subtraction method 

Of the different tested parameters we find that the stamp selec- 
tion method affects the subtraction quality and detection success 
the most. The automatic stamp selection in ISIS mostly works 
fine, but the lack of control in selecting stamps can be problem- 
atic when dealing with real SN survey data, due to the presence 
of saturated objects and the rejection routine sometimes failing 
for varying objects. Using a set of sub-optimal stamps will result 
in a subtracted frame with worse subtraction quality. We note 
that the optimal selection for stamps is based on SNR, selecting 
preferentially point-like sources as stamps, i.e. in general fainter 
sources than in the SNR based selection, degrades the subtrac- 
tion quality. Varying the widths of the Gaussian basis functions 
of the kernel within a reasonable range does not significantly al- 
ter the detection results, however it should be noted that if all 
of the functions are too wide, depending on the data in ques- 
tion the image subtraction will fail. We also note that for surveys 
with small seeing differences, like the SVISS, the default ker- 
nel widths used in ISIS might be to wide. Using widths more in 
line with the theoretically optimal convolution kernel width can 
marginally improve the subtraction quality. 



6.4. Implications for SN surveys 

The detection efficiency of supernovae in a multi-epoch sur- 
vey can be studied by Monte Carlo simulations of SN light 
curves as observed by the telescope used for the survey (see e.g. 
iDahlen et al. 1 12004 . The light curves, together with the individ- 
ual detection efficiencies for each epoch can then be used to find 
an overall detection efficiency for supernovae of a certain type 
and at a given redshift. The flux offset effect found in this work 
will affect the detection limits as well as the number of spuri- 
ous detections. These can be handled in different ways, but their 
presence will always be a problem and minimising the number of 
spurious detections is frequently important. A positive flux offset 
will result in an artificially boosted detection efficiency and an 
increase in the number of spurious detections. A negative offset 
will instead cause the detection efficiency to go down, making 
the survey to lose depth. This effect can possibly be countered 
by controlling the detection threshold used by SE using a weight 
map, in principle redefining the threshold for positions consid- 
ered to be inside galaxies. 

It is very unlikely that systematic flux offsets of the kind 
found in this work affect the results found by the large 
SN la projects. In these projects (e.g . iBarris & Tonrvl 120061: 
lAstier et al.]|2006t iMiknaitis etaLll2007L flie photometric accu- 
racy is studied extensively to make sure that the systematic er- 
rors in the cosmological parameter determination is minimised. 
As shown above the detection efficiency will be influenced by 
the flux offset. Projects studying the supernova rate and doing 
detections in images subtracted using methods similar to ours 
might suffer from this effect and detailed simulations should be 
run to make sure that no offset is present (or correct for it other- 
wise). The results from some recent surveys ( Dahlen et al.l2(j04t 
ICappellaro et"!!] 120051: IBarris & Tonrvll2006l) are based on su- 
pernovae with observed spectroscopy. The limiting magnitudes 
for these surveys are thus set by the magnitude limit of the spec- 
troscopic instrument used , which is considerably brighter than 
the corresponding imaging limit, thus the offset problem de- 
scribed here is likely not an issue. However, with present and 
upcoming surveys based on imaging alone and pushing the de- 
tections to the fainter limits this effect needs to be understood 
and studied thoroughly. 



Based on the findings reported in this paper we can make 
some recommendations for supernova detection and photome- 
try using image subtraction methods. Testing should be done on 
at least part of the data to determine optimal parameters for the 
image subtraction. The stamp selection should be controlled in 
some way, a good option is to create a master stamp list from 
the reference image and use that in all subtractions but making 
sure that none of the stamps contain SNe. If the automatic stamp 
selection in ISIS is used, care must be taken that the selected 
stamps are not variable or contain saturated pixels. To make sure 
that a possible flux offset will not affect the results and to es- 
timate photometric errors, simulations should be done. If pho- 
tometry is done in the subtracted frames, the measurements for 
an individual SN should be carefully done. This can be done by 
simulating SNe of similar magnitude in nearby galaxies of the 
same brightness and size as the host galaxy and by measuring 
the possible offset. The photometric errors for each point on the 
light curve should be estimated by using the magnitude scatter 
found from simulations. If no offset correction for the detection 
efficiencies is applied we recommend that the detection efficien- 
cies should be capped by the efficiency for isolated SNe (in the 
presence of a positive flux offset), the true efficiency can never 
be better inside than outside galaxies. 



7. Summary 

We have described the SN detection and photometry pipeline 
used for the SVISS. We have also presented detailed testing of 
supernova detections and photometry using image subtraction 
and automatic detection. Our main results and conclusion are: 

- Image subtraction using ISIS 2.2 works well for the SVISS, 
the subtraction quality as measured by noise levels in the 
location of galaxies is in general good. Simulations show that 
our pipeline is successful in detecting SNe as well as in doing 
photometry on them. 

- The detection efficiency and photometric accuracy for su- 
pernovae in a single epoch is affected by how the stamp se- 
lection for ISIS is done. We have also found that the host 
galaxy brightness does affect these quantities, likely through 
the presence of powerful subtraction residuals at the location 
of bright galaxies. 

- Image subtraction using ISIS benefits from using a con- 
trolled stamp list. We have found that the best subtraction 
quality is achieved when the subtraction stamps used have 
high signal-to-noise ratios, ideally the brightest non-varying 
sources in the field should be used as stamps. The number of 
stamps depends on the size and quality of the image, and for 
our data set ~60 stamps per (Ikxlk) pixels has proven to be 
sufficient. 

- We have discovered a systematic flux offset in subtracted im- 
ages. This offset only appears at the locations of galaxies and 
is quite small. Photometry of simulated SNe in isolated po- 
sitions show no offset. For individual sources the offset is 
smaller than the estimated photometric scatter. This effect 
can cause systematic errors in the photometry if not cor- 
rected for Detection limits of SNe will also be affected by 
non-zero offset, causing problems with increased numbers 
of spurious detections for positive offsets and a decrease in 
survey depth for a negative offset. It should be noted that, 
based on the tests run so far, we cannot say for sure whether 
the offset is something unique to the SVISS data set or if it 
is due to the image subtraction method. 
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- Simulations of supernovae both in and outside galaxies are 
needed to obtain proper estimates of the photometric errors. 
The simulations can also be used to correct for possible sys- 
tematic offsets. 

- Using automatic source detection (in this case SE) we find 
that the detection efficiencies can be pushed to the limits of 
the data. The number of spurious detections in a single epoch 
subtraction is very large, on the order of 80%, but decreases 
as a function of subtraction quality. To reject spurious detec- 
tions other methods must be used (e.g. multiple filters, light 
curve fitting). 

Understanding possible systematic effects of the subtraction and 
detection process is very important when interpreting the results 
from the SVISS, but will also be imperative for the upcoming 
supernova projects based on large imaging surveys, e.g. Pan- 
STARR^i LSSlEI and SNAP (lAlderindboOSl) . 
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